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ABSTRACT 



We deal with a Yukawa-like long-range modified model of gravity (MOG) 
which recently allowed to successfully accommodate many astrophysical and cos- 
mological features without resorting to dark matter. On Solar System scales 
MOG predicts retrograde secular precessions of the planetary longitudes of the 
perihelia w whose existence has been put on the test here by taking the ratios 
of the observationally estimated Pitjeva's corrections to the standard Newto- 
nian/Einsteinian perihelion precessions for different pairs of planets. It turns 
out that MOG, in the present form which turned out to be phenomenologically 
successful on astrophysical scales, is ruled out at more than 3o" level in the Solar 
System. If and when other teams of astronomers will independently estimate 
their own corrections to the usual precessions of the perihelia it will be possible 
to repeat such a test. 

Subject headings: Experimental tests of gravitational theories; Modified theories of 
gravity; Celestial mechanics; Orbit determination and improvement; Ephemerides, 
almanacs, and calendars 
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Introduction 



The modified gravity (MOG) tlieory put forth in fjMoffatll2006l ) was used successfully to 



describe various observational 
resorting to dark matter (see 



phenomena on a s 



troph ysical and cosmological scales without 



Moffat and TothI (120081 ) and references therein). It is a fully 



covariant theory of gravity which is based on the existence of a massive vector field coupled 
universally to matter. The theory yields a Yukawa-like modification of gravity with three 
constants which, in the most general case, are running; they are present in the theory's 
action as scalar fields which represent the gravitational constant, the vector field coupling 
constant and the vector field mass. Actually, the issue of the running of the parameters 



of modified models of gravity is an old one, 



inown in somewhat similar contexts since the 



early 1990s (see, e.g., (IBertolami et aL 



1993f) and reference s therein). An approximate 



20071 ) allows to compute their values 



solution of the MOG field equations (IMoffat and Toth 
as functions of the source's mass. 

The resulting Yukawa-type modificat ion of the inverse-square New ton's law in the 



gravitational field of a central mass M is f lMoffat and Toth 



2007 



20081) 



A 



MOG 



{1 + a [1 — (1 + /ir) exp(— /ir)]} 



(1) 



where Gn is the Newtonian gravitational constant and (IMoffat and Toth 



2007, 



20081 ) 



M 



a 



M + C[ 



G 



- 1 , Goo ^ 20 Gn, Gi ^ 25000 M^ 



G2 ^ 6250 Mq/^ kpc"\ 



Such values have been obtained by (IMoffat and Toth 



1/2 

' 



(2) 



(3) 



20071 ) as a result of the fit of the 



velocity rotation curves of some galaxies in the framework of the searches for an explanation 
of the fiat rotation curves of galaxies without resorting to dark matter. 
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In this paper we will put eq. ([T]) on the test in the Solar System in order to check 
if it is compatible with the latest determinations of the corrections (A ri?) to the usual 



Newtonian/Einsteinian planetary perihelion precessions (jPitjeva 



2005a 



b|). Such quantities 



are solved-for parameters phenomenologically estimated in a least-square sense by fitting 
about one century of observations of various kinds concerning planetary motions with their 
values theoretically predicted by simultaneously integrating the equations of motions of 
the major bodies of the Solar System written in terms of a complete suite of dynamical 
force models. They encompass all the known Newtonian e ffects along with those pr edicted 



by general relativity, except for the Lense-Thirring effect (ILense and Thirring 



19181 ) which 



was left unmodelled (see below). Thus, (Arij) account for any unmodelled/mismodelled 
dynamical features of motion and are suitable, in principle, to put on the test MOG, 
provided that they are used in a suitable way. 



Note that Moffat and TothI (120081 ) explicitly write that eq. ([T]), with eq. and eq. 
([3]), is not in contradiction with the present-day knowledge of Solar System dynamics. We 
will show that it is not so also for any other (non-zero) values of a and /i, with the only quite 
general condition that /ir <^ 1 in Solar System, as it must be for any long-range modified 
model of gravity. It is interesting to point out that Yukawa -like modifications of Newton' s 



20051). 



law might also be tested in the context of stellar dynamics (IBertolami and Paramos 
Here we outline the procedure that we will follow. 

Generally speaking, let LRMOG (Long-Range Modified Model of Gravity) be a given 
exotic model of modified gravity parameterized in terms of, say, K, in a such a way that 
K = would imply no modifications of gravity at all. Let P(LRMOG) be the prediction of 
a certain effect induced by such a model hke, e.g., the secular precession of the perihelion 



- 5 - 



of a planet. For all the exotic models considered it turns out thac 

V (LRMOG) = Kg{a,e), (4) 

where g is a. function of the system's orbital parameters a (semimajor axis) and e 
(eccentricity); such g is a. peculiar consequence of the model LRMOG (and of all 
other models of its class with the same spatial variability). Now, let us take the ratio 
of 'P(LRMOG) for two different systems A and B, e.g. two Solar System's planets: 
Pa(LRMOG)/Pb (LRMOG) = gx/QB- The model's parameter K has now been canceled, 
but we still have a prediction that retains a peculiar signature of that model, i.e. g^l Qb- Of 
course, such a prediction is valid if we assume K is not zero, which is just the case both 
theoretically (LRMOG is such that should K be zero, no modifications of gravity at all 
occurred) and observationally because K is usually determined by other independent long- 
range astrophysical/cosmological observations. Otherwise, one would have the meaningless 
prediction 0/0. The case K = {] [oi K < K) can be, instead, usually tested by taking one 
perihelion precession at a time. If we have observational determinations O for A and B of 
the effect considered above such that they are affected alscG by LRMOG (it is just the case 
for the purely phenomenologically estimated corrections to the standard Newton-Einstein 
perihelion precessions, since LRMOG has not been included in the dynamical force models 
of the ephemerides adjusted to the planetary data in the least-square parameters' estimation 



process by Pitjeva (jPitjeva 



2005a 



|b|)), we can construct Oa/C^b and compare it with the 
prediction for it by LRMOG, i.e. with gh/gB- Note that 50/0 > 1 only means that 
O is compatible with zero, being possible a nonzero value smaller than 60. Thus, it is 



^In our case it will he K = — a/i^, as we will see in Section [2J 

^If they are differential quantities constructed by contrasting observations to predictions 
obtained by analytical force models of canonical Newtonian/Einsteinian effects, O are, in 
principle, affected also by the mismodelling in them. 
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perfectly meaningful to construct Oa/O-q- Its uncertainty will be conservatively evaluated 
as \1/0b\50x + \0x/0'^\50-Q. As a result, Ca/C^b will be compatible with zero. Now, the 
question is: Is it the same for g^/ Qb as well? If yes, i.e. if 

Ca Pa(LRMOG) 



Cb Pb(LRMOG) 
within the errors, or, equivalently, if 

Ca Pa(LRMOG) 



(5) 



(6) 



Cb Pb(LRMOG) 

within the errors, LRMOG survives (and the use of the single perihelion precessions can be 
used to put upper bounds on K). Otherwise, LRMOG is ruled out. 

2. The predicted perihelion precessions and the confrontation with the 

measured non-standard rates 

In the case of the Sun, eq. (j2]) and eq. (j3]) yield 

00 ^ 3 X 10-^ /i ^ 3 X 10"^ AU-\ (7) 

so that 

ttQ/i^ = 3 X 10~^^ AU~l (8) 

Since in the Solar System /ir 10~^ — 10~^, we can safely assume exp(— /ir) ^ 1 — /ir, so 
that 

Amog^-^(1 + «/xV). (9) 
As a result, a radial, uniform perturbing acceleration 

A = -G^Mafx^ ^ 10"^^ m s'^ (10) 

is induced. 
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The secular, i.e. averaged over one orbital revolution, effect of a small radial and 
unform perturbing accelerat i on on the longitude of the perihelion of a planet w has been 
worked out by, e.g., ISandersI (120061 ): it amounts to 



dw \ . ail — e^l 



irh^\/^-"^Wa(l-e^), (11) 

Clearly, using only one perihelion rate at a time would yield no useful information on MOG 
due to the extreme smallness of the perturbing acceleration, as told us by eq. ( ITOl) . Thus, let 
us take the ratios of the perihelion precessions. It must be noted that the following analysis 
is, in fact, truly independent of the values of a and /i, provided only that afi'^r'^ <^ 1 in the 
Solar System so as that the perturbative approach can be applied to eq. ([9]); the condition 
/ir <^ 1 is the cornerstone of any long-range modified models of gravity, and should a ^ 1 
the planetary orbits would have been distorted in a so huge manner that it would have been 
detected since long time. Applying the scheme outlined in Section [T] to our case in which 
K = — a/i^ and g{a, e) = \/ G^iMaiX — e^), one can construct 

n = ^ / 12 

with the estimated corrections {^zu) to the standard Newtonian/Einsteinian perihelion 
precessions of planets A and B, listed in Table [T|, and compare them to the theoretical 
prediction 

A^J^^^. (13) 

y flBll - e|) 

obtained from eq. (11 ip . for that pair of planets A and B. The results are in Table [2J |n — ^| 
is different from zero at more than 3cr level for A = Venus, B = Mercury, A = Earth, B = 
Mercury and A = Mars, B = Mercury. It is important to note that the errors have been 
conservatively evaluated as 

< ini f + i<^) (14) 
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Table 1: Inner planets. First row: estimated corrections to the standard precessions of the 



longi tudes o 



3 of (jPitjeva 



the pe rihelia in 10 ' cy (" cy — > arcseconds per century), from Table 



2005bl ) (apart from Venus). The quoted errors, in 10~^ " cy~^, are not the 
formal ones but are realistic. The formal errors are quoted in square brackets (E.V. Pitjeva, 
personal communication to the author, November 2005). The units are 10~^ " cy~^. Second 
row: semimajo r axes, in Astronomical Units (AU). Their formal errors are in Table IV of 



(jPitjeva 



2005al ). in m. Third row: eccentricities. Fourth row: orbital periods in years. The 
result for Venus have been recently obtained by including the Magellan radiometric data 
(E.V. Pitjeva, personal communication to the author, June 2008). 





Mercury 


Venus 


Earth 


Mars 


{Aw) (10-4 " cy-i) 


-36 ± 50 [42] 


-4±5[1] 


-2 ±4[1] 


1±5[1] 


a (AU) 


0.387 


0.723 


1.000 


1.523 


e 


0.2056 


0.0067 


0.0167 


0.0934 


P (yr) 


0.24 


0.61 


1.00 


1.88 



Table 2: First column: pair of planets. Second column: 11 for that pair of planets. The 
errors come from the realistic uncertainties in (Azu). Third column: A for that pair of 
planets. Fourth column: a level of discrepancy between 11 and A for that pair of planets. 



A B n A (X 

Venus Mercury 0.1 ±0.3 1.4 4 

Earth Mercury 0.05 ±0.18 1.64 8 

Mars Mercury -0.03 ±0.18 2.02 11 



J 



because of the existing correlationqfl among the estimated corrections to the precessions of 



^The maximum correlation, 26%, occurs for the Earth and Mercury (E.V. Pitjeva, per- 
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perihelia. 



If we repeat our analysis by subtracting from (Azu) the main canonical unmodelle d 



effect, i.e. the general relativistic Lense-T 
induced by the Sun's angular momentum (llorio 



lirring prec essions ( iLense and Thirring 



19181) 



2007al ) shown in Table [3l i.e. if we use 



(15) 



(A^b)* {Awb) - 4'^^ ' 
the situation does not substantially change, apart from the sigma level at which |n* — ^| is 
not compatible with zero, as shown in Table HI 

Table 3: General relativistic Lense-Thirring precessions of the longitudes of perihelia of the 
inner planets of the Solar System in 10~^ " cy~^. 



Mercury 


Venus 


Earth Mars 


^(LT) (10-4 " cy-1) -20 


-3 


-1 -0.3 



Table 4: First column: pair of planets. Second column: 11* for that pair of planets including 
the unmodelled general relativistic Lense-Thirring effect. The errors come from the realistic 
uncertainties in {Avo) . Third column: A for that pair of planets. Fourth column: a level of 
discrepancy between 11 and A for that pair of planets. 

A B n* A a 

Venus Mercury 0.06 ± 0.51 1.4 2.7 
Earth Mercury 0.06 ± 0.44 1.64 3.5 
Mars Mercury -0.08 ± 0.56 2.02 3.7 



sonal communication to the author, November 2005). 
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The availability of the corrections to the usual rates of perihelia of several planets 



allows us to put on the test MOG also in another way as well. The acceleration 



([T]) can also be recast in the commonly used Yukawa form (IMoffat and Toth 

GyM 



20071 ) 



aw of eq. 



A 



Y 



r\ / r' 
1 + aY ( 1 + exp y-jj^ , 



where 



G 



G, 



Y 



ay 



1 + ay 
(Goo - Gn)M 



(a 



Gn)A/ + Gn( 
1 



M + c[y' 



A 



In the case of the Sun 



a 







-3.04 X 10-^ Gy = 1.00000003040Gn, A = 33000 AU. 



(16) 

(17) 
(18) 
(19) 

(20) 



loriol fl2007br ) in 



A Yukawa-type acceleration of the form of eq. (fT6!) has been tested by 
the Solar System without any a-priori assumption on the size o^ ay; concerning A, it was 
only assumed that A > ae. By using the correctio ns to t he stan dard rates of the perihelia 



of A = Earth and B = Mercury quoted in Table [T] 



loriol fl2007bf ) found 



A 



In 



OA 



0.182 ±0.183 AU, 



(21) 



which contradicts eq. (1201) . Using the data for Venus in the equation for Q;v (jlorio 



2007bh 



yield; 



2A2 (Aw) 
ay = — , exp , , 

VGyMa VA 



(22) 



ay = (-1 ±4) X 10 



-11 



(23) 



"^The strength parameter a used in (jloridl2007bl ) can be identified with ay here. 

^According to eq. fl20l) . using Gn in eq. (122|) instead of Gy, as done in (llorid l2007bl ). 
does not produce appreciable modifications of the results. 
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which is three orders of magnitude smaller than the result of eq. (120|) . 



If we use n* for the Earth and Mercury in eq. fl2T]) and (Azu)* for Venus in eq. 
the results does not change appreciably; indeed, we have 

A = 0.2 ± 0.4 AU, ay = (-0.3 ± 2.7) x 10^^^ 



(24) 



3. Conclusions 



In the framework of the attempts of explai ning certain a strophysical and cosmo logical 



2006 



Moffat and Toth 



2003) is 



features without invoking dark matter, MOG (iMoffat 
a long-range modified model of gravity, based on a vector field and three scalar fields 
representing running constants, which assu mes a Yukawa-like form . Rece nt developments of 



this theory allowed their proponents to fix flMoffat and Toth 



2007 



20081 ) the values of the 



constants entering it. We have shown that, on Solar System scales, MOG yields a uniform 
anomalous acceleration which would induce retrograde planetary perihelion precessions. 
We put on the test the possibility that such exotic precessions exist by comparing the 
ratio of them A for different pairs of planets to the ratio 11 of the corrections to the 
usual Newtonian/Einsteinian precessions estimated by E.V. Pitjeva which account for any 
unmodelled/mismodeleld dynamical effects. It turns out that 11 7^ ^ at more than 3a level 
even by including in 11 the main unmodelled canonical effect, i.e. the general relativistic 
Lense-Thirring precessions. Conversely, using the estimated corrections to the planetary 
perihelion rates to phenomenologically determine the strength parameter of the putative 
MOG Yu kawa force and its range yie. 



of MOG (iMoffat and Toth 



2007 



2008 



ds values which are neatly incompatible with those 
B|). In assessing the results presented here it must be 
considered that, at present, no other people have estimated the non-standard part of the 
planetary perihelion motions; it would certainly be useful to repeat the present analysis 



if and when other teams of astronomers will estimate their own set of corrections to the 
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standard perihelion precessions as well. 
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